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^ ■ Using the Generalized Maximum Entropy Principle based on the 

o ■ 

! nonextensive q entropy a new family of random matrix ensembles is gen- 

erated. This family unifies previous extensions of Random Matrix Theory 
(RMT) and gives rise to an orthogonal invariant stable Levy ensemble 
with new statistical properties. Some of them are analytically derived. 
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O ■ Random matrix theory (RMT) started in physics with the introduction by E. 



Wigner, in the 50s, of Gaussian matrix ensembles, the Orthogonal (GOE), the Unitary 



■ GUE) and the Symplectic (GSE). Their properties were fully developed by Dyson, 
>• . Gaudin, Mehta and others [1]. These ensembles have a wide application as models to 

X 

H ' describe statistical properties of quantum fluctuations of systems of few or many-body 

particles. They have been useful in discussing nuclear and atomic properties, meso- 
scopic physics, quantum chaos, theory of amorphous solids, etc (see, for instance [2]). 
The link between RMT and Information Theory was set by Balian [3] who, by using 
the Boltzman-Gibbs-Shannon entropy associated to the ensemble probability distribu- 
tion, obtained the Wigner ensembles by maximizing it subjected to the normalization 
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condition and a constraint given by the average norm of the matrices. Ensembles to 
describe symmetry breaking have been constructed by adding an extra constraint to 
this scheme [4]. 

In this letter, we use this framework and consider ensembles within the General- 
ized Maximum Entropy Principle (GMEP) based on the nonextensive Tsallis entropy 
[5]. This entropy has been applied to a great variety of phenomena, specially those 
in which long range correlations are present (see however [6] concerning its physical 
interpretation). It is dependent on the non-additivity parameter q defined in such a 
way that when q ^ 1 the Boltzman-Gibbs-Shannon entropy is recovered. We show 
that a new family of ensembles is generated that unifies some important extensions of 
RMT. In the range — oo < ^ < 1, it is found to be a restricted trace ensemble that 
interpolates between the bounded trace ensemble [7] when g — > — oo and the Wigner- 
Gaussian ensembles &X, q — 1. In the domain 1 < q < q'max, with gfmax being a cutoff 
imposed by the normalization condition, it interpolates between RMT at g = 1 and an 
ensemble of Levy matrices [8] that appears at the neighborhood of the extremum g^ax 
where the ensemble distribution has divergent moments. 

As extensions of RMT that preserve the stability of the universal ensembles. Levy 
matrices have attracted recently much attention due to its potential application to 
many areas ranging from physics to finances [8-10]. Stability means that if Hi and 
H2 are matrices of the ensemble, their sum H — Hi + H2 also is [11]. This will be 
the case if the individual matrix elements are distributed according to a Gaussian or 
a Levy function. We prove that this indeed happens, in the case of the g-generalized 
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ensembles, for all allowed values of q, i.e. — oo < q < qmax when N goes to infinity. 

Although the individual matrix element distribution of g-ensembles have the same 
asymptotic behavior as the Levy matrices of Ref . [8] , there is here a basic difference 
as they are orthogonal invariant with matrix elements, in principle, correlated. Or- 
thogonal invariance is also satisfied by the ensembles of Refs. [9,10] which are directly 
defined in terms of the joint distribution of eigenvalues. However no explicit reference 
to the matrix elements distribution is made there and the spectral statistical measures 
are obtained expressing them in terms of apropriately defined orthogonal polynomials. 
Here we do not apply this technique and show that the special relation that g-ensembles 
have with the Gaussian ensembles allows their spectral properties to be analytically 
derived. 

Applied to matrices whose entries are random variables, the nonextensive entropy 
can be written as 



where H is a, NxN matrix distributed according to P (H) and dH is the product of dif- 
ferentials of the independent variables of the matrices. For definiteness we consider real 
symmetric matrices in which case we have / = A'^(A^-|-l)/2 independent matrix elements 




1 - / dHP" {H) 



(1) 



and the differential in (1) is conveniently defined as dH — 2^^^ ^^''^ ni<i<j<Ar 



dHi^. 



The GMEP consists in maximizing (1) subjected to normalization 




(2) 



and to the constraint [12] 
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j dHP'i (H) tTH^-fxj dHPi (H) = (3) 

that fixes the gf-average of the norm defined as the trace of the square of the matri- 
ces. Following the usual steps of the variational method, we arrive at the probability 
distribution 

P{H;\a)^Z],'(l + jtrH''^~' (4) 

with A given by 

X — — au — r — — . 5) 

q-1 ^ q-1 2 ^ ^ 

Zjv in (4) is the partition function and (3) has been used to determine the relation 
Oi — ^ ■ Let us remark that had we used Renyi's entropy [13] instead of Eq. (1) we 
would also have been led to Eq. (4). 

Changing from matrix elements to eigenvalue and eigenvector variables the ensemble 
distribution factorizes and, after integrating over the eigenvector parameters, we find 
for the eigenvalues the joint probability distribution 

P{E,, ...E^; A, a) = Xjv (l + ^ E ^l) ~ U \Ej - E^l , (6) 

where Kj^ is the normalization constant. Taking — > 1 in the above, A — > oo and the 
RMT distributions 

Pgoe (H; a) = Zq^^ j^ exp (-a tr if^) (7) 

and 
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Pgoe{Eu -En; a) = Kgoe,n exp (-a ^ UlEj - Ei 



(8) 



are recovered. 

Considering — oo < g < 1, i.e. — | > A > — oo the condition tr(i7^) — Y^El< — ^ 
has to be imposed in order to warrant a real positive probabihty distribution for any 
q. These two inequahties define hyperspheres in which the matrix elements and the 
eigenvalues are confined in their respective spaces. Taking in Eq. (4) the limit q —>■ — oo 
with the partition function given by 



where © (x) is the step function. The bounded trace ensemble is known to follow the 
Wigner- Dyson statistics of the Gaussian ensemble when — > oo. To show that this 
is also the case for — oo < g < 1 we consider the probability distribution of a generic 
matrix element 




(9) 



we find that the ensemble goes to the bounded trace ensemble 




(10) 





and the correlation between two matrix elements hi and /i2 




(12) 



By taking the hmit of large matrices, (11) goes to the Gaussian distribution 
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p {x; A, a) ~ \ — exp (— ax^) (13) 

V TT ^ ^ 

while C (/ii, /12) — > indicating that the matrix elements behave as those of the Gaus- 
sian ensembles as A?" — > 00. Numerical simulations [14] confirm that the level density 
is given by the Wigner semi- circle law 



Pgoe (E; a) 



2^ N_^2 \E\<J^ 
■K \J a ' II \J a 



(14) 



0, \E\ > 
and spectral fluctuations follow GOE statistics. 

Consider now q> 1. The partition function is given by 

that requires the restriction A > or g < gmax = 1 + 1- We see that the introduction of 
the parameter A is crucial to be able to study the limit N ^ 00. It maps the interval 
1 < q < Quiax onto the interval 00 > A > 0. The Fourier transform of the distribution 
of a generic matrix element, Eq. (11), with A > is 

A-.(*y'f) (16) 

where Kx{z) is the modified Bessel function [15]. In order to ensure that spectra 
scale independently of the size of the matrices, a has to go to infinity when N does. 
This can be seen from the analytic expression of the level density, Eq. (27) below. 
The requirement is that a characteristic value, say E^. = \J^, remains finite when 
diverges. In this limit, Kx (z) can be replaced by its small z expansion and keeping 



only the first terms we can write F {k; A, a) ~ exp (^~^ 



k /A 

2 V a 



with 



(7 = 2 and A 



1 

4(A^ 



if oo > A > 1. 



(17) 



and 



(7 = 2A and A = ^ — ^ if 1 > A > 

r(i + A) 



(18) 



Therefore for oo > A > 1 the distribution of a generic matrix element approaches the 
Gaussian distribution 



, X , /(A-l)a 
p {x; A, a) ~ y — — exp 



a 



- (A - 1) jx' 



(19) 



For 1 > A > the Levy-Gnedenko generahzed central hmit holds [16] and p {x; A, a) 
goes to the Levy function, L^. {x, a, A) = tt"^ /q°° (it exp {—Af^) cos {xt), with the same 
asymptotic behavior, i.e. 



a 



p{x;X,a) ~ 2;/-L2A 



' fa r(l-A) 



(20) 



Concerning correlations between matrix elements, Eq. (12) shows that only for large 
values of A or (3; the matrix elements behave independently, whereas for small values, 
A < 2, they are strongly correlated. Therefore for large X or a (19) predicts that the 
level density goes to the semi-circle paoE E; (A — 1) ^ i.e. Eq. (14) with a replaced 
by(A-l)f. 

We focus now on the spectral properties of these new ensembles. They are analyt- 
ically derived by introducing the representation 



a 



A 
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, (ieexp(-Oe^"'exp(-^etri/' 



(21) 



that allows the joint distribution function of the matrix elements to be written in terms 
of the joint distribution function of the GOE ensemble as 

P {H- A, a) = ^ di exp (-^ C'-'^goe [h- j ; (22) 

the joint distribution of eigenvalues becomes 
P (E,, E^; X, a) = exp (-0 C^"' exp (-^^ E ^l) U \Ej - E,\ 

(23) 

where Xjv is the normalization constant. Integrating (23) over all eigenvalues we deduce 
the relation 

Kn = r(A) "^^ ^^^'^ ^^^^ 

relating Kn to the corresponding RMT constant in standard units, i.e. a = | in Eq. 
(8) see [1]. Substituting in (23) one finally obtains for the normalized joint eigenvalue 
density 

N_ 

P {E^, ...E^- A, «) = ^ ) ' d^exp (-0 e^'^-'PaoE (^^i, ...^x^r; ^) 



(25) 



where we have introduced the rescaled eigenvalues — y^E^. This is one of the 
central results of this paper and can be taken as the defining equation of the new 
ensemble. It expresses the eigenvalue distribution of the new ensemble as a sort of 
r function of the GOE eigenvalue distribution. It shows that one may expect that 
measures of the g'-family will be weighted Laplace transforms of the corresponding 
measures of the Gaussian ensemble. 
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Integrating (25) over all eigenvalues but one and multiplying by N, the average 
eigenvalue density is expressed in terms of Wigner's semi-circle law as 



(26) 



The asymptotic power law behavior of this distribution is better seen by rewriting it 



as 



p{E;\,a) 



N 



\E 



2A+1 



f— I .T7TW^^fA+2,A + 2,--^ (27) 



J r(A)r(A + 2) 



where M (a, b, z) is the confluent hypergeometric function [15]. In Fig. 1, with a — ^ 
(see Eqs. (17) and (18)) the density p {E; A, a) is plotted for four values of A, exhibiting 
the deviation from the semi-circle law as A moves inside the interval 1 > A > 0. When 
A — > 0, the density behaves as p ~ ,Jl\+i approaching the same behavior as for a 
nonconfining log square potential [17]. 

The behavior of the spectral fluctuations can be illustrated by considering the gap 
probability function E (s) (usually denoted E (0, s)) that gives the probabihty of finding 
an eigenvalue-free segment of length s. This function has been investigated in Ref. [10] 
for Cauchy ensembles and is related to the presence of gaps in the spectrum. For the 
g-family it is expressed in terms of the corresponding GOE function as 



1 r'^ 



A 
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(28) 



obtained integrating the joint eigenvalue density over all eigenvalues outside the interval 
(—6*, 9) around the origin. In (28) y {x) = 2 /q dtpaoE {t) . Together with 



s(e)^2 f dEp (E; A, a) 
Jo 



(29) 



(28) expresses E (s) in a parametric form. Using the Wigner surmise for the nearest 
neighbor spacing distribution p (s) and the relation connecting E (s) and p (s), Eqoe 
in (28) can be well approximated by Egoe (?/) — !— erf On Fig. 2 results in 

the Levy regime are displayed. Notice the large increase of the probability of formation 
of a gap with respect to the GOE case. The asymptotic behavior in Eq. (28) can be 
extracted by making the substitution x — \f^0 that leads to 



2a \9 



x^^-'EGOE[y{x)]. (30) 



For large 6*, this equation predicts for A = 1 a power law decay E (s) ~ clearly seen 
in the figure. This very characteristic behavior is exhibitted here for the first time. 

In summary, we have proved that the g-generalized family of ensembles interpolates 
between the bounded trace ensemble [7] at the extremum — > — 00 and the Wigner- 
Gaussian ensembles at = 1. In the domain 1 < q < gmax, it interpolates between 
RMT at g = 1 and an ensemble of Levy matrices at the neighborhood of the extremum 
?max = 1 + J- These orthogonal invariant stable matrix ensembles have novel spectral 
properties. Remarkably, several of their distribution functions can be expressed as 
integral transforms (sort of extended F functions) of the corresponding distribution 
functions of the Gaussian ensembles. 

It is premature to exhibit specific applications of these generalized ensembles. How- 
ever there are worth exploring possibilities, for instance, connections with the so-called 
critical statistics [18] or the transition from Erdos-Renyi to scale free models in random 
graph theory [19]. In conclusion, let us remind that stable laws (Levy laws) were first 
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introduced and studied. It was correctly anticipated [16] that a large domain of appli- 
cations would follow [20]. We believe that we are presently facing a similar situation, 
where the role of a random variable is now being extended to the one of a random 
matrix. The results presented here should contribute to broaden the applications of 
random matrix theory. 

After completion of this letter, we learned of ref. [21] closely related to the work 
presented here. 

Fruitful discussions with C. TsaUis are acknowledged. A.C.B. and M.P.P. are sup- 
ported by the Conselho Nacional de Pesquisas (CNPq). This work is supported by a 
project CAPES-COFECUB. 
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Figure Captions 

Fig. 1 The eigenvalue density for four values of the parameter A (= 10, 1, 0.75, 0.5) 
in the transition region from the Gaussian to the Levy regime, with N — 50. For the 
sake of comparison the semi-circle Pgoe E; (A — 1) ^ 
shown (dashed hue). 

Fig. 2 The eigenvalue-free probability E (s) for A = 1. Full line: theory, Eq. (28) 
and its asymptotics (dotted line); dashed line: Egoe '■ numerical simulation with 
N — 20. See text for further explanation. 
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with A = 10 and a = ^ is also 
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